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Abstract 

Intracellular recordings of cortical neurons in vivo display intense subthreshold membrane po- 
tential (V„,) activity. The power spectral density (PSD) of the V„, displays a power-law structure 
at high frequencies (>50 Hz) with a slope of about -2.5. This type of frequency scaling cannot 
be accounted for by traditional models, as either single-compartment models or models based on 
reconstructed cell morphologies display a frequency scaling with a slope close to -4. This slope is 
due to the fact that the membrane resistance is "short-circuited" by the capacitance for high fre- 
quencies, a situation which may not be realistic. Here, we integrate non-ideal capacitors in cable 
equations to reflect the fact that the capacitance cannot be charged instantaneously. We show that 
the resulting "non-ideal" cable model can be solved analytically using Fourier transforms. Nu- 
merical simulations using a ball-and-stick model yield membrane potential activity with similar 
frequency scaling as in the experiments. We also discuss the consequences of using non-ideal ca- 
pacitors on other cellular properties such as the transmission of high frequencies, which is boosted 
in non-ideal cables, or voltage attenuation in dendrites. These results suggest that cable equations 
based on non-ideal capacitors should be used to capture the behavior of neuronal membranes at 
high frequencies. 
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1 Introduction 



One of the greatest achievements of computational neuroscience has been the development of cable 
theory (reviewed in [[Il|2l), and which can explain many of the passive properties of neurons, including 
how dendritic events are filtered by the cable structure of dendrites. Cable theory describes the space 
and time propagation of the membrane potential by partial differential equations. Such a formalism 
constitutes the basis of nearly all of today's computational models of dendrites, and is simulated by 
several publically-available and widely-used simulation environments (reviewed in [|3l). 

Some experimental observations, however, may suggest that the standard cable formalism may not be 
adequate to simulate the fine details of dendritic filtering. One of these observations is the fact that 
the power spectral density (PSD) of synaptic background activity or channel noise does not match 
that predicted from cable theory |l4l[5l[6l|7]|. The PSD scales approximately as 1/f" with an exponent 
a = 2.5, both for channel noise and background activity (Fig. [iK-B), whereas cable theory would 
predict scaling with an exponent a = 4 or a = 5 for synaptic inputs distributed in dendrites [l5l[8]|; see 
also Appendix 1), or a = 3.2 to 3.4 when inputs are distributed in soma and dendrites (see Fig. \TiZ- 
D). In other words, these data suggest that frequencies are filtered by dendritic structures in a way 
different from that predicted by traditional cable equations. 

Figure [H here 

One possible origin of such a mismatch could be due to the fact that the permittivity of the membrane 
is frequency dependent [HKTOl. However, capacitance measurements in bilipid membranes shows neg- 
ligible variations around 100 Hz (see Fig. 5 in ifTOl ). suggesting that the Cole-Cole model may not be 
the correct explanation for this range of frequencies. It could also be that distortions of the frequency 
dependence arise from the complex three-dimensional morphology of the neuronal membrane [fTTI . 
However, NEURON simulations of the standard cable model using three-dimensional morphologies 
of cortical pyramidal neurons give frequency scaling with an exponent a > 3 (Fig.[T]Z!-D), suggesting 
that this is not a satisfactory explanation either. 

None of the previous models take into account the fact that the surface of neuronal membranes is a 
complex arrangement, not only of phospholipids, but also of a wide diversity of surface molecules [fT2l|. 
This complex surface may be responsible for additional resistive phenomena not taken into account 
in previous approaches. In other words, the neuronal membrane may not be an "ideal" capacitor, as 
commonly assumed in the standard cable formalism. In the present paper, we explore this hypoth- 
esis as an alternative mechanism to explain the observed frequency scaling and consider neuronal 
membranes as "non-ideal" capacitors. We show that cable equations can be extended by including a 
non-ideal resistive component (Maxwell- Wagner time) in the capacitor representing the membrane, 
and that the non-ideal cable model reproduces the observed frequency scaling. We also show con- 
sequences of this extension to cable equations in voltage attenuation and synaptic summation. Our 
aim is to provide an extended cable formalism which is more adapted to capture membrane poten- 



tial dynamics and dendritic filtering at high frequencies. Some of these results have appeared in a 
conference abstract [[T3l . 



2 Materials and Methods 



The standard and non-ideal cable equations were either solved analytically (see Results) or simulated 
using custom-made programs written in MATLAB. A "ball-and-stick" model consisting of a soma 
connected to a dendritic cylinder of length 1^ was simulated (see Results for details). Away from the 
current source, we have the following equations (in Fourier space): 

X'^^^^^ = KL(a))V,„(x,a)) (1) 

2 ' ^ - 



1 + /COXm 



where X = \fr^^i is the electrotonic constant that characterizes the cable, x,„ is the membrane time 
constant, and Xm is the Maxwell-Wagner time constant (Xm = corresponds to the standard cable 
equations; see Results). 

The "source" synaptic current consisted in a random synaptic bombardment of Poisson-distributed 
synaptic events. Each synaptic event consisted of an instantaneously rising current followed by expo- 
nential decay, and were summated linearly: 

h= AY^Hit- 1,) exp[-(? - t,)lxs\ , (2) 

i 

where /5 stands for the source current, H{t) is the Heaviside function, and ti are the times of each 
synaptic event (Poisson-distributed with mean rate of 100 Hz). The decay time constant was X5 = 
10 ms and the amplitude of the current was A = 1 nA. 

The source current was inserted at different positions in the dendrite (see Results). The voltage at the 
soma was obtained by solving either standard or non-ideal versions of cable equations (see Results 
and Appendix 2). The power spectral density (PSD) was calculated from the somatic membrane 
potential using the fast Fourier transform algorithms present in MATLAB (Signal Analysis toolbox). 
The same algorithm was also used to calculate the PSD from experimental data. 

The experimental PSD of V,„ activity shown here were obtained from intracellular recordings of cat 
parietal cortex neurons in vivo and were taken from previous publications [|4ll2l|, where all method- 
ological details were given. No filter was used during digitization of the data, except for a low-pass 
filter with 5 kHz cutoff frequency during acquisition (sampling frequency of 10 kHz). Thus, the PSD 
is expected to reflect the real power spectral content of recorded V,,, up to frequencies of 4-5 kHz. 



Some simulations (Fig. [i|C-D) were realized using morphologically reconstructed neurons from cat 



cortex obtained from two previous studies [fT4l [TSl . where all biological details were given. The 
three-dimensional morphology of the reconstructed neurons was incorporated into the NEURON 
simulation environment, which enables the simulation of the traditional cable equations using a three- 
dimensional structure with a controlled level of spatial accuracy lfT6l . Simulations of up to 3500 
compartments were used. In vzvo-like activity was simulated using a previously published model of 
synaptic bombardment at excitatory and inhibitory synapses [il7il (see this paper fo details about the 
numerical simulations). 

3 Results 

We start by deriving the non-ideal cable model, then investigate its general properties by evaluating 
the PSD of somatic voltage, as well as voltage attenuation. 

3.1 Derivation of non-ideal cable equations 
3.1.1 The membrane as a non-ideal capacitor 

In electrostatics, if an electric field is applied to a closed conductive surface, electric charges mi- 
grate until they reach equilibrium (when the field tangential to the surface is zero). In particular, 
the electric resistivity of the membrane imposes a given velocity to charge movement, which dissi- 
pates calorific energy similar to a friction phenomenon. This calorific dissipation is usually neglected, 
which amounts to consider an instantaneous charge re- arrangement following changes in electric field. 

However, in reality this calorific dissipation may have significant consequences, and this phenomenon 
is well known for capacitors |18| . A "non-ideal" capacitor dissipates calorific energy when the elec- 
tric potential varies, and capacitors are usually conceived such as to minimize this phenomenon and 
realize the well-known ideal relation i = C^. A "non-ideal" linear capacitor can be represented as an 
arrangement of resistances, inductance and capacitance (see Fig.[2K)- A linear approximation, which 
is usually sufficient for most purposes. In particular, this approximation is valid when the effects of 
electrostriction are negligible [[T0l[T9l . This is the case when the propagated signals are of small am- 
plitude (millivolts), because C{V) = C(0) (1 +aV^), with typically a = 0.02 V^^ [IgJ- In such cases, 
the membrane capacitance can be represented by a resistance and a capacitance in series ll20l (see 
Fig.[2&). The resistance represents here the loss of calorific energy associated with charge movement. 
In standard cable equations, such a resistance is not present (see Fig.[2lC). 

Thus, we use a more realistic capacitor modeled by taking into account an additional resistance (Rsc), 
which accounts for the calorific loss and the consequent finite-velocity of charge rearrangement. This 
R-C circuit will be characterized by a relaxation time Xm = RscC, called "Maxwell time" or "Maxwell- 
Wagner time" II2TII22I . The Maxwell time corresponds to the characteristic displacement time of the 



charges in the capacitor. Thus, such a non-ideal capacitor cannot be charged instantaneously; the 
resistance Rsc imposes a minimal charging time due to finite charge velocities. 



Figure [2] here 

This phenomenon of finite charge velocity is particularly relevant to biological membranes, which are 
capacitors in which charges are also subject to rearrangements. In the following, we attempt to include 
this contribution to membrane capacitors by including Maxwell- Wagner time to cable equations and 
determine its consequences. 



3.1.2 Non-ideal cable equations 

We extend cable equations by including a finite charge velocity (or equivalently, a minimal charging 
time) to membrane capacitors. We start by Ohm's law, according to which the axial current ii in a 
cylindric cable can be written as: 

i-^oE^ ^ . (3) 
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We also have, for the membrane current i 
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and we can write 

.„=^+r^^w (5) 
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where Cm(?) is the inverse complex Fourier transform of the capacitance c,„((0). c,n{t) = c,n5(f) if the 
capacitance does not depend on the frequency. 

Integrating Maxwell- Wagner phenomena, we have: 
Thus, we obtain the following non-ideal cable equations: 



1^ "^""^'■"'y- dt ^^-^^^^^ 

dCn,{t-t') 



where X = \fh^i is the electrotonic constant that characterizes the cable. 



3.1.3 General solution of non-ideal cable equations 



The non-ideal cable equations (Eqs. [6l) are a linear system with constant coefficients which can be 
solved by using Complex Fourier Transforms: 

/oo 
V,n{x,t) e'"^' dt 
-oo 

/oo 
Vc{x,t) e'"" dt 
-oo 

/oo 
Cn,{t) e'"' dt 
-oo 

We obtain the following expression: 

" kLv„(x, (0) (7) 



with 
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where Xm((o) = r,nC,„i(i>) and Xm((o) = r5cC,„((o) are the parameters that characterize the cable. 
The general solution of Eq. |7]is given by 

v,„(;c,a)) = A(a))exp( ^-"^^ "-'^ )+B{oy)exp{- ^-^^^ ) (9) 
where Is is the position of the current source in the dendrite. 

This solution is similar to that of traditional cable equation, with the only difference in the value of K. 
In cable equations, this value is given by 

= 1 + /COXm . (10) 

In particular, for null frequency, the two cable formalisms are equivalent 

K..f(0)=K,(0) = l , (11) 

whereas they will predict different behavior for (O > 0. 

In the following, we will consider that the capacitance is independent of frequency, Cm{(o) = est, as 
also assumed in the standard cable model E El- 
Figure [3] compares the values of K between the two cable formalisms (with Cm(co) = est). The differ- 
ence depends on the relative values of Xm and x,„: for Xm « x,„, the two formalisms are very similar, 
but differ when Xm is larger, in particular for high frequencies. Thus, the critical parameter is Xm, 
which determines the saturation of the value of K. 



Figure [3] here 



3.2 Voltage attenuation vs. distance and frequency 

To compare the properties of the non-ideal cable model compared to the standard cable model, we 
evaluated the properties of voltage attenuation in a large dendritic branch. We have chosen a cable of 
Id = 500 /jm and diameter of 2 /jm, with a current source situated at one end of the cable (x = Is = 0) 
and connected to an infinite impedance at the other end (x = Z^; "sealed end"). In these conditions, 
we can determine the law of voltage attenuation with distance, using complex Fourier analysis. 

As we have seen above, the main difference between the standard and non-ideal cable models lies in 
the expression for K (see Eqs.[8]and[T0l). In a finite cable of constant diameter, the steady-state voltage 
attenuation profile is given by the relation: 

Vi„{x,0)) = A((o) exp{--x) +B{(o) exp{-x). (12) 

k A 

for X > 0. To evaluate the functions A(ol)) and 5((o), we apply the limit conditions of the dendrite. 
At X = 0, we have a current source is = I = id, and at x = Id we have id = ("sealed end"). The 
expressions for A and B are then given by Eqs.[19]and[20l respectively (see Appendix 2). 

This relation is plotted in Fig.|4]for two values of the membrane time constant x,„ of 5 ms and 20 ms, 
which correspond to two different conductance states of the membrane (the corresponding electro- 
tonic constant is A = 353.5 /jm and 707.1 /jm, respectively). The voltage attenuation is in general 
steeper for the non-ideal cable model, which effect is particularly apparent for frequencies of the or- 
der of 0-50 Hz. However, this effect reverses between 50 and 100 Hz, in which case the non-ideal 
cable model shows a less steep voltage attenuation profile compared to the standard cable model (see 
50 and 100 Hz in Fig.®. 

Figure |4] here 

3.3 Power spectra of voltage noise predicted by non-ideal cable equations 

We now calculate the PSD of the voltage noise predicted by non-ideal cable equations. We consider 
a "ball-and-stick" model consisting of a soma and a dendritic segment of variable length (Fig. [5K)- 
The source consists of a sum of exponentially-decaying currents (see Materials and Methods), which 
represent the synaptic current resulting from many synapses releasing randomly, as shown in Fig. [5^. 
The source has a PSD which scales as 1 //" with an exponent a = 2 at high frequencies (Fig.[5]C!). 



Figure [5] here 



To investigate the PSD of the somatic voltage in the ball-and-stick model, we first examine the PSD 
following a single source consisting of summated exponential synaptic currents. The standard cable 
model predicts that such a source localized on a dendritic branch (ball-and-stick model with 1^ = 
500 /jm and X ~ 400 ^m) gives a V,„ PSD scaling approximately as 1//" with an exponent a ~ 
4, which corresponds to a somatic impedance much larger than that of the dendrite (soma radius 
of 7.5 /jm; see Appendix 1), which would correspond to most central neurons for which the soma 
represents a minor proportion of the membrane. The V„, PSD for the standard cable model with 
uniformly distributed exponential synaptic currents is illustrated in Fig. [6] (continuous curve), and 
shows a frequency scaling with an exponent a ~ 4. 

In contrast, the non-ideal cable model gives different scaling properties of the PSD, according to the 
value of Xm (Fig. [6l dotted and dashed lines). The power for high frequencies (>50 Hz) is much larger 
in the non-ideal cable model compared to the standard model, which shows that non-ideal cables have 
enhanced signal propagation for high frequencies. The V„, PSD for the non-ideal cable model with 
uniformly distributed exponential synaptic currents is illustrated in Fig. [6] (dashed curve), and shows 
a frequency scaling with an exponent 2 < a < 4 for x,„ > Xm > 0, respectively (a ~ 2 when x,„ = Xm, 
but it can be shown that a = 2 only if Xm ^ °°)- 

Figure [6] here 

We next investigated the influence of the localization of the current source in the dendrite. Figure [TK 
shows the PSD obtained at the soma of the ball-and-stick model when the current source was placed at 
different positions in the dendrite. The position affects the amplitude of the PSD, and the frequency- 
scaling of the PSD is affected by the position. The scaling exponents obtained are of a = 4.1416 for 
250 fjm and 5.3653 for 450 ^m for the standard model, and a = 2.5311 for 250 ^m and 2.8354 for 
450 /jm for the non-ideal cable model. The PSD obtained when simulating a "distributed" synaptic 
bombardment in the dendrite (Figure |7^) also displays the same frequency- scaling. Similar results 
were also obtained by varying the parameters x„j and Xm (not shown), suggesting that the properties 
of frequency scaling, as shown in Fig.[6l are generic. 

Figure |7] here 

To evaluate the optimal value of Xm (for this particular model with Xm = 5 ms), we fitted the PSD of the 
model to that of experiments. To perform this fit, we used a frequency range of 100 to 400 Hz, which 
was chosen such that it is not affected by instrumental noise (<700 Hz) and such that the frequency 
band considered belongs to the power-law scaling region of the spectra (>80 Hz). The result of this 
fitting is shown in Fig. [H The scaling exponent obtained are of a = 3.6533 for the standard cable 
model, and of a = 2.3306 for the non-ideal cable model, for an optimal value of Xm = 0.3 x,„. This 
suggests that the calorific dissipation caused by the resistivity of the membrane to charge movement 



is of the order of 30% of that caused by the flow of ions through ion channels. This estimate is of 
course specific to the model used, but variations of this model {Id, diameter, number of dendrites, for 
a uniform Xm over the whole neuronal surface) showed little variation around this value (not shown). 

Figure [8] here 

This value gives a cutoff frequency (1/ Xm) around 105 Hz. Above this cutoff frequency, the membrane 
becomes more resistive than capacitive because the energy loss due to calorific dissipation becomes 
larger than the energy necessary for charge displacement. This is very different than an ideal capacitor, 
in which the energy from the current source would exclusively serve to charge displacement. In Fig.[3l 
one can see that the value of k for the non-ideal model departs from that of the standard cable model 
around this cutoff frequency. 

Thus, from the above figures, and especially Fig. |6l it is apparent that the non-ideal cable model has 
more transmitted power compared to the standard cable model at high frequencies (>>100 Hz). This 
increased transmission of high frequencies is also visible by superimposing the V,,, activities of the 
standard and non-ideal model (Fig. |9l). Such an increased transmission at high frequencies can be 
explained by the fact that in the standard cable model, the term 1 / z(Oc„j tends to zero when CO tends 
to infinity, such that for high frequencies r„, is short-circuited by the capacitance of the membrane. 
In the non-ideal cable model, such a short-circuit does not occur, even at frequencies much larger 
than the cut-off frequency. This results in a very different behavior at high frequencies, and a less 
pronounced frequency fall-off in the non-ideal cable PSD. Displacing charges by capacitive effect 
takes energy, and this energy diminishes with increasing frequencies in the non-ideal cable, which 
enables more energy transfer between remote ion channels in dendrites (synapses for example) and 
the soma at high frequencies. This is also consistent with the fact that the non-ideal cable equations 
display less voltage attenuation (see Section IX2l) . 

Figure [9] here 

4 Discussion 

In the present paper, we have proposed an extension to the classic cable theory to account for the 
behavior of neuronal membranes at high frequencies. Experimental observations indicate that the 
PSD of the V„, does not match that predicted from cable theory, in particular for the frequency- 
scaling at high frequencies [|4l|5l|6ll3. The modification to cable equations consists in incorporating 
a "non-ideal" membrane capacitance by taking into account the calorific dissipation due to charge 
displacement, which is usually neglected. We have shown that this "non-ideal" cable formalism can 
account for the frequency scaling of the PSD observed experimentally for high frequencies (Fig. [8]). 



In experiments with channel noise or synaptic noise, the PSD scales as 1/ with an exponent a 
around 2.5 [ffl [51 [6l |3 . The standard cable model predicts that the somatic V,,, should scale with an 
exponent a comprised between 3 and 4 |l5l, when the source is located in the soma. However, we 
have shown here that the frequency scaling of the V,„ PSD depends on the location of the source, and 
that the exponent a is equal or larger when current sources are located in dendrites (see Fig. |7] and 
Appendix 1). Thus, the standard cable model cannot account for exponents lower than a = 3. On the 
other hand, taking into account non-ideal capacitances may lead to scaling exponents down to a = 2, 
depending on the magnitude of the dissipation in the non-ideal capacitance (as quantified by the value 
of the Maxwell- Wagner time Xm', see Fig. [6]). In the case Xm is non-uniform, then one may have larger 
differences of frequency scaling between somatic and dendritic current sources (not shown). 

In the non-ideal model, the calorific dissipation originates mostly from the resistance of the mem- 
brane to lateral ion displacement. This "tangential" resistance is not yet characterized experimentally 
and is equivalent to the resistance involved in the non-instantaneous character of membrane polar- 
ization Il22]| . Several arguments indicate that this resistance may be substantial. First, the membrane 
surface contains various molecules such as sugars and various macromolecules, in addition to phos- 
pholipids [[T2|. Thus, lateral ion movement is likely to be affected by collisions or tortuosity imposed 
by these molecules. Second, the phospholipids themselves contain local dipoles at their polar end, 
which is likely to cause local electrostatic interactions which may influence the lateral movement of 
ions. Indeed, the fitting to experimental data using the non-ideal cable model predicts a value for Xm 
which is a significant fraction (~30%) of the membrane time constant. 

The complex three-dimensional membrane morphology could have consequences on frequency-de- 
pendent properties even with traditional cable theory ifTTI . We tested this possibility by simulating 
detailed three-dimensional morphological models of cortical pyramidal neurons and failed to repro- 
duce the frequency scaling of the activity in vivo (see Fig. [U). Thus, although the morphology 
does affect frequency scaling, it does not account for the values observed experimentally. 

Another source of distortion in the frequency dependence of the V,„ is the fact that membrane per- 
mittivity (and capacitance) may also depend on frequency [|9l|23l. Such a frequency dependence is 
caused by a calorific dissipation during the polarization of the membrane flU, while the Maxwell- 
Wagner phenomenon that we discuss here is a calorific dissipation during the movement of charges 
on the membrane surface. However, direct capacitance measurements of bilipid membranes do not 
evidence any significant variation of permittivity for frequencies around 100 Hz ifTOl . and thus cannot 
explain the observed deviations between cable theory and experiments shown in Fig. \T\ Moreover, 
these measurements ifTOl [T9l were realized on artificially reconstructed membranes, which have a 
much simpler structure compared to neuronal membranes (no saccharides, no proteins, etc). This 
is compatible with the possibility that in biological membranes, the Maxwell- Wagner effect may be 
particularly prominent. The dependence of the membrane capacitance c„, on frequency may explain 
the flattening of the PSD above 1000 Hz, which is visible in the experimental PSDs (see Fig. [8]). How- 
ever, the most likely explanation for this flattening is that the recording is dominated by instrumental 
noise at such frequencies (note that the bending of the experimental PSD above 4000 Hz in Fig. [8] is 



likely due to the low-pass 5 kHz filter used during data acquisition). 

Other factors may also affect the frequency scaling. Taking into account the finite rise time of synap- 
tic events by using double exponential templates amounts to add a factor 2 to the exponent a [[HI. 
Similarly, introducing correlations in the presynaptic activity may also affect the frequency scaling 
of V„, power spectra |[24l . Li all these cases, however, the change in the scaling always consists in 
increasing the exponent a, while a decrease is needed to account for a=2.5 scaling. 

Thus, the frequency scaling of the V,,, activity can be affected by several factors as discussed above. 
Our results show that the non-ideal character of the neuronal membrane can account for the observed 
frequency scaling. We believe that in reality, a combination of factors is responsible for the observed 
frequency scaling, and future experiments should be designed to test which are the most determinant 
on frequency scaling, and what are the consequences on the integrative properties of neuronal cable 
structures. 

Finally, our results show that the frequency-dependence of the steady-state voltage profile (Fig. [4]) is 
also affected by the non-ideal character of the membrane capacitance. Simulations show that high- 
frequency signals (> 100 Hz) propagate over larger distances in the non-ideal cable model compared 
to the standard cable model. This theoretical result may be important to understand the propagation 
of high-frequency events such as the "ripples" oscillations [|25ll26l across dendritic structures. 

In conclusion, we provided here an extension to cable equations which incorporates the non-ideal 
character of the membrane capacitance. We showed that this extension yields several detectable con- 
sequences on neurons. First, it affects basic cable properties such as the voltage attenuation profile, es- 
pecially at high frequencies. Second, it radically changes the frequency- scaling properties of voltage 
power spectra. The observed frequency scaling is within the range predicted by the non-ideal cable 
model. Fitting the model to experiments provides an estimate of how "non-ideal" is the membrane 
capacitance, and the significant values of Xm found here suggest that indeed, neuronal membranes 
may be far from being ideal capacitors. 



Appendix 1: Frequency scaling in the standard cable model 



In this Appendix, we overview the frequency scaling characteristics of the PSD of the V,,, for the 
ball-and-stick model using the standard cable equations. 

Dendritic current source located close to the soma 

We first consider the ball-and-stick model with an isolated current source located in the dendrite close 
to the soma. From expression (l24l) (see Appendix 2), we have: 

(Z2 ©Z3)/^o ~ lim(Z2 ©Z3)) = Z3 , 

and from expression (12T]) . when the distance / from the source to the soma is small, the impedance of 
the distal part of the dendrite is given by 

Zi ^ —coth{^)] , 

where Id is the length of the dendrite. From expression (fT4l) . for small /, we have 

VE = FAis ^ (— II Z3)/5 , 

where ^ is the input impedance of a finite dendritic branch. Thus, from expression (1281) . for small /, 
we obtain 

Fr(Z,(o) ^limFr(Z,(0) = 1 . 
Because Fg ~ 1, the membrane potential at the center of the soma is given by 

Vsoma = ( II Zi,)is (13) 

when the current source is located close to the soma. 

Thus, for high frequencies (> 100 Hz), the PSD of the somatic V,,, scales as 1/ with a G ]3,4[ for 
a exponential current source located close to the soma. This result is similar to single-compartment 
models d. 

General case of dendritic current source 

We now consider the general case of a current source located at an arbitrary position in the den- 
dritic branch of the ball-and-stick model. We have necessarily Ft ^ I, resulting in a supplementary 
dependence on frequency. Moreover, the current divider Fa also depends on frequency. Numerical 



simulations show that the PSD of the somatic V,„ scales as 1/ with an exponent a > 3. For ex- 
ample, with exponential currents uniformly distributed on a dendrite of 1^ = 500 fjm, the frequency 
scaling is close to an exponent of a=4 (see continuous curve in Fig. (6]). We verified numerically (not 
shown) that the standard cable model cannot give a frequency scaling with a slope smaller than a = 3 
(using Poisson-distributed synaptic inputs). 

A similar scaling with an exponent a = 4 was observed earlier, when simulating realistic dendritic 
morphologies based on reconstructed cortical pyramidal neurons 



Appendix 2: Impedance analysis of the ball-and-stick model 

In this appendix, we derive the expressions needed to study the frequency dependence of the ball- 
and-stick model (Fig.[5K), for both standard and non-ideal cable equations. The ball-and-stick model 
consists of a soma, which is assumed to be the recording site, and a dendritic branch which contains 
the source. Referring to Fig. [5K, we have the source (S) and the recording locations (P), as well as 
the impedances corresponding to the different regions (Zi for the distal part of the dendrite, away of 
the source, Z2 for the proximal part of the dendrite, between the source and the soma, and Z3 for the 
soma). 

We first evaluate the voltage at the current source: 

. Zi(Z2ffiZ3) . 

'^■^'' Z,+(Z.ffiZ3) ='^'"- 

where the term (Z2 ©Z3) is the input impedance of the dendritic segment in series with Z3. Fa is 
the input impedance as seen by the current source is located at a position Is on the dendritic branch. 
Expression [141 shows how Fa varies as a function of the position of the source in the dendrite. 

Next, we calculate the somatic voltage from the transfer function of the dendritic branch. Ft, which 
links the voltage at the source with the somatic voltage. 

Vso,na=FTVE (15) 



Finally, we calculate the voltage transferred to the soma from the equivalent circuit (Fig.fTOl). 

Vp = — ^^^;=—Vsoma = Fb Vsoma 1 (16) 

where Fb is the voltage divider caused by the fact that the tip of the recording pipette is located inside 
the soma at some distance from the membrane (in case of sharp-electrode recordings). This divider is 
entirely resistive and very close to 1, which expresses the fact that the exact position of the pipette is 
not a determining factor in the value of Vp. 



Figure [To! here 

Thus, we have 

^soma = Fb Ft Fa is ^ Ft Fa is (17) 
We calculate these different terms below. 

Input impedance Zi (distal part of the dendrite) 

For a current source is located at position Is, we have 

is = idi {h, w) + id2 (h, (O) , (18) 

where ij^ (Z^, Co) is the current density at the beginning of the distal part of the dendrite (of length 
All) , and id2{ls,(ii) is the current density of the proximal part of the dendrite (see Fig(fTOl). From 
expression [9l we have 

r; ax kri 



where |£| > can be as small as desired. This factor arises because we consider point current sources, 
in which case the spatial derivative of the Vm is discontinuous at ;c = /,,. 

From the "sealed end" condition, we have 

u(4 + A/i,a)) = --^(/^ + A/i,a)) = [A((o) exp{-^^^) -5((o) exp(^^)] = . 

r,- ax KYi K K 



Thus, we have 



and 



Consequently, we obtain 



l-expi^) 



Zi = ^-7z T = — coth[^—) (21) 



where k = or 'Kg^t for standard or non-ideal cable models. 



Input impedance (Z2 0Z3) (proximal region) 



For the proximal part of the dendrite (of length A/2 = h), which is in series with the impedance Z3 at 
x — (see Fig. [To]), we have (see expressions [TSl and |9]) 

. n \ 1 3v,„(/,- |£|,a)) K 

where |e| > can be as small as desired. 
Moreover, we have 



and 



Thus, we obtain 



and 



Consequently, we obtain 



ri dx Xri X A, Z3 



kZ kZ 
v„,(0,(o) = A((o) exp(^) +5((0) 



5(a))=A(a)) 



5((o) =A((o) + — Zrf2(^.^'«) 



Afro) = Xr/ [KZ3->ir,]Zrf,(Z„a)) 

K [exp(^) + 1][kZ3 ?an/z(f ) +;ir,-] 

and 

Xr,[KZ3 + Xr,-]Z^,(Z,,(D)e.yp(^) 
K [exp(^) + 1][kZ3 z'an/z(f^) +Xr,-] 

Thus, the input impedance (Z2 ©Z3) is given by: 

(Z2©Z3) ^ = ^^^4^ + X^r^tanhj^^) 

id, (Z„ CO) [KZ3 tank In] K [KZ3 tank ( ^ ) + ?ir,] 

where X = and k = or Ke;^; according to which cable model is used. 
For Z3 00, we obtain the input impedance from Eq. [ 



Calculation of the transfer function Fj 



To evaluate Ft, we calculate the voltage at points = / by imposing Vmijs, (o) = 1 at points = Is. With 
this initial value, the voltage v,„(x) at point = equals the value of the transfer function at point 
x — Q (see Eq. (9]). In such conditions, we obtain: 



A(a))+5((o) = 1 . 



Thus, we have 

IC K K 

Fr(x,(0) =A((0) [exp{-{ls~x))- exp{--{ls-x))]+exp{--{ls-x)) (25) 
The voltage v„, at point x = must equal Z3 /,(/, (o) (current conservation). We have 

-r; li 



dx 

Consequently, we must obtain 



^r— U=o= -ri — x=0= Tl V,n x=0= Tl x=0 (26) 
ox Z3 

where ^1 = — 27- Thus, we have 

^ {^-M)exp{-'f) 

^[exp{^^) + exp{-'^)]+Xx\ [exp{f)-expi-f)] 

and the transfer function is given by 

Fr(0,a)) =A(a)) [exp(|/,) - exp{-^h)]+exp{-^ls) (28) 



Finally, we have 

Z3 = Z30 + . — , p N , 1 (29) 

where 230 is the plasma resistance in the soma. K equals Kv or Kgxt according to the cable model 
considered. 



References 



[1] Rail, W. 1995. The Theoretical Foundation of Dendritic Function (Segev, I., J. Rinzel and G.M. Shepherd, 
ed). MIT Press, Cambridge MA. 

[2] Johnston, D. and S.M. Wu. 1995. Foundations of Cellular Neurophysiology, MIT Press, Cambridge MA. 

[3] Brette, R., M. Rudolph, T. Carnevale, M. Hines, D. Beeman, J.M. Bower, M. Diesmann, A. Monison, 
P.H. Goodman, F.C. HaiTis Jr., M. Zirpe, T. Natschlager, D. Pecevski, B. Ermentrout, M. Djurfeldt, A. 
Lansner, O. Rochel, T. Vieville, E. MuUer, A. Davison, S. El Boustani and A. Destexhe. 2007. Simulation 
of networks of spiking neurons: A review of tools and strategies. J. Computational NeuroscL, in press. 



(article available at http : //arxiv . org/abs/q-bio . NC/06110891 



[4] Destexhe A., M. Rudolph and D. Pare. 2003. The high-conductance state of neocortical neurons in vivo. 
Nature Reviews Neurosci. 4: 739-751. 

[5] Diba, K., H.A. Lester and C. Koch. 2004. Intrinsic noise in cultured hippocampal neurons: experiment 
and modeling. /. Neurosci. 24: 9723-9733. 

[6] Jacobson, G.A., K. Diba, A. Yaron-Jakoubovitch, C. Koch, I. Segev I and Y. Yarom. 2005. Subthreshold 
voltage noise of rat neocortical pyramidal neurones. /. Physiol. 564: 145-160. 

[7] Rudolph M., J.G. Pelletier, D. Pare and A. Destexhe. 2005. Characterization of synaptic conductances and 
integrative properties during electrically-induced EEG-activated states in neocortical neurons in vivo. /. 
Neurophysiol. 94: 2805-2821. 

[8] Destexhe, A. and M. Rudolph. 2004. Extracting information from the power spectrum of synaptic noise. 
/. Computational Neurosci. 17: 327-345. 

[9] Cole, K.S. and R.H. Cole. 1941. Dispersion and absorption in dielectrics. I. Alternating cunent charac- 
teristics. /. Chem. Phys. 9: 341-351. 

[10] White, S.N. 1970. A study of lipid bilayer membrane stability using precise measurements of specific 
capacitance. Biophys. J. 10: 1127-1148. 

[1 1] Eisenberg, R.S. and R.T Mathias. 1980. Structural analysis of electrical properties. Crit. Reviews Bioeng. 
4: 203-232. 

[12] Alberts, B, A Johnson, J Lewis, M Raff, K Roberts and P Walter. 2002. Molecular Biology of the Cell, 
Fourth Edition. Garland Publishing, New York. 

[13] Destexhe, A. and C. Bedard. 2007. A non-ideal cable formalism which accounts for fractional power- 
law frequency scaling of membrane potential activity of cortical neurons. Soc. Neurosci. Abstracts 33: 
251.14. 

[14] Contreras, D., A. Destexhe, and M. Steriade. 1997. Intracellular and computational characterization of 
the intracortical inhibitory control of synchronized thalamic inputs in vivo. J. Neurophysiol. 78: 335-350. 



[15] Douglas R.J., K.A. Martin and D. Whitteridge. 1991. An intracellular analysis of the visual responses of 
neurones in cat visual cortex. /. Physiol. 440: 659-696. 



[16] Hines, M.L. and N.T. Carnevale. 1997. The NEURON simulation environment. Neural Computation 9: 
1179-1209. 

[17] Destexhe, A. and D. Pare. 1999. Impact of network activity on the integrative properties of neocortical 
pyramidal neurons in vivo. /. Neurophysiol. 81: 1531-1547. 

[18] Bowick, C. 1982. RF Circuit Design, Newnes Elsevier, New York. 

[19] Alvarez O. and R. LatoiTC. 1978. t Voltage-dependent capacitance in lipid bilayers made from monolayers. 
Biophys. J. 21: 1-17. 

[20] Raghuram, R. 1990. Computer Simulation of Electronic Circuits, John Wiley & Sons, New York. 
[21] Raju, G.G. 2003. Dielectrics in Electric Fields, CRC Press, New York. 

[22] Bedard, C, A. Destexhe and H. Kroger. 2006. Model of low -pass filtering of local field potentials in brain 
tissue. Physical Review E73: 05 19 1 1 . 

[23] Hanai T., D.A. Haydon and J. Taylor. 1965. Some further experiments on bimolecular lipid membranes. 
/. Gen. Physiol. 48 (suppl): 59-63. 

[24] Marre, O., S. El Boustani, P. Baudot, M. Levy, C. Monier, N. Huguet, M. Pananceau, J. Fournier, A. 
Destexhe and Y. Fregnac. 2007. Stimulus-dependency of spectral scaling laws in VI synaptic activity as 
a read-out of the effective network topology. Soc. Neurosci. Abstracts 33: 790.6. 

[25] Ylinen, A., A. Bragin, Z. Nadasdy, G. Jando, I. Szabo, A. Sik and G. Buzsaki. 1995. Sharp wave- 
associated high-frequency oscillation (200 Hz) in the intact hippocampus: network and intracellular 
mechanisms. /. Neurosci. 15: 30-46. 

[26] Grenier, F, I. Timofeev and M. Steriade. 2001. Focal synchronization of ripples (80-200 Hz) in neocortex 
and their neuronal correlates. /. Neurophysiol 86: 1884-1898. 



Figure Legends 



Figure 1: Fall-off structure of power spectra of synaptic noise in cortical neurons. A. Time course 
of the membrane potential during electrically-induced active states in a cortical neuron recorded in- 
tracellularly from cat parietal cortex in vivo (data from [Q)- B. Power spectral density (PSD) of the 
membrane potential in log scale. The PSD has a fall-off structure which follows a power law with 
a fractional exponent, around -2.6 in this case (dashed line; modified from refs. |4l|7]|)- C. Four dif- 
ferent morphologies of cortical pyramidal neurons from cats obtained from previous studies [fT4l[T5l . 
and which were incorporated into numerical simulations. D. PSD obtained from the four models in 
C, using the traditional cable formalism in NEURON simulations. The power-law exponent obtained 
was of 3.4, 3.3, 3.2 and 3.4, respectively (cells shown from left to right in C). 



Figure 2: Different equivalent electric schemes for capacitors. A. Linear model of a capacitor, con- 
sisting of two resistances (Rsc and Rpc), one inductance (L^c) and one capacitance element (C). B. Ap- 
proximation of the linear model obtained by including a resistance (Rsc) in series with the capacitance 
(C). This leads to a characteristic relaxation time for charging the capacitor (given by Xm = RscC). C. 
Ideal capacitance as in the standard cable model. 



Figure 3: Comparison between K values in the standard and non-ideal cable model. The values of 
K are plotted for the two models for various values of Xm and two values of x,,, (5 ms and 20 ms). 
The function K saturates for the non-ideal cable model, and the value of the saturation equals to 
y^l +t„t/xm- The K curves for the non-ideal model depart from the standard model for a frequency 
that approaches the cut-off frequency of /c = 1/ 2j^- 



Figure 4: Steady-state voltage profile in a finite cable. A cable of 500 /urn length and 2 /urn diameter 
was considered with a current source atx = {Cm = 1 /jF/cm^; R,- = 2 Q.m). The voltage profiles in the 
non-ideal (gray lines) and standard (black lines) cable models are compared for different frequencies. 
Two values of the membrane time constant are considered, x,,, = 5 ms (A) and x^ = 20 ms (B), which 
correspond to two different conductance states (Xm = 1 .5 ms in both cases, which corresponds to Xm 
= 0.3 x„i in A, and Xm = 0.075 x,„ in B). 



Figure 5: Ball-and-stick model used for calculations. A. Scheme of the ball-and-stick model where 
P indicates the soma, S the position of the current source, and Z1...Z3 are impedances used in the 
calculation. B. Example of a source current representing synaptic bombardment in the ball-and-stick 
model. The current source consists in Poisson-distributed exponential currents (see Materials and 
Methods). C. Power spectral density of the synaptic current source shown in B. The PSD scales as a 
Lorentzian (1//" with an exponent a = 2 between 100 and 400 Hz). 



Figure 6: Power spectral density of the V,„ of the ball-and-stick model with exponential synaptic 
currents uniformly distributed in the dendrite (from 1 to 450 /j m, every 10 /j m). The current source of 
each synaptic event was the same and equals exp(— f /0. 1 ) nA, and the PSD is shown for the membrane 
potential at the soma. The continuous curve shows the standard cable model, while the other curves 
(dotted and dashed) show the non-ideal cable model with different values of Tm- Parameter values: 
C,„ = 1 A/F/cm^, x„, = 5 ms. Id = 500 /jm, = 1 ^um, Rsoma = 7.5 /jm, R, = 2 Q.m. 



Figure 7: Power spectral density of multiple synaptic events in the ball-and-stick model. A. Voltage 
PSD at the some for a source current similar to Fig.[5B which was placed at different positions in the 
dendrite (from top to bottom: 250 and 450 ^um from the soma). For each location, the PSD is shown 
for the standard cable model (gray) and for the non-ideal cable model (black). B. PSD obtained when 
the source currents were distributed in the dendrite (from 1 to 450 /urn, every 10 /urn). Parameter 
values: C,„ = 1 /jF/cm^, x„, = 5 ms, = 500 /urn, Rj = 1 /jm, Rsoma = 7.5 /jm, R,- = 2 Q.m, Xm = 0.3 x„i- 



Figure 8: Best fit of the non-ideal cable model to the power spectral density obtained from intracel- 
lular experiments. The non-ideal cable model was simulated using a ball-and-stick model subject to 
synaptic bombardment (see Materials and Methods). The dendritic branch had a 75 /Jin length and 
the power spectral density (PSD) was calculated from the somatic membrane potential. Black: ex- 
perimental PSD (see Fig.[T]); Gray: model PSD (see Fig.[5]C! for the PSD of the current source). The 
slopes were calculated using a linear regression in the frequency band 100-400 Hz. The optimal value 
for Xm was of 0.3x,„. Parameter values: C,„ = 1 /jF/cm^, x,,, = 5 ms, = 75 /urn, Rd = I A^m, Rsoma = 
7.5 /jm, Rj = 2 Q.m. 



Figure 9: Comparison of V„, activities in the standard and non-ideal cable models. The current source 
is indicated on top, while the bottom trace shows the V,„ activities superimposed. The inset shows a 
detail at 5 times higher temporal resolution. Same parameters as the optimal fit in Fig.[8l 



Figure 10: Equivalent circuit for the ball-and-stick model. Zi is the input impedance of the dendritic 
branch (open circuit), Z2 is the impedance of the intermediate segment, in series with the impedance 
Z3 of the soma. 
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